Magnetization reversal via internal spin waves in magnetic nanoparticles 
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By numerically solving the equations of motion for atomic spins we show that internal spin- 
wave processes in large enough magnetic particles, initially in unstable states, lead to complete 
magnetization reversal and thermalization. The particle's magnetization strongly decreases in the 
middle of reversal and then recovers. We identify two main scenarios, exponential and linear spin- 
wave instabilities. For the latter, the longitudinal and transverse relaxation rates have been obtained 
analytically. Orientation dependence of these rates leads to a nonexponential relaxation of the 
particle's magnetization at long times. 

PACS numbers: 75.50.Tt, 75.75.-|-a, 75.30.Ds, 76.20.-|-q 



I. INTRODUCTION 

From both fundamental and application viewpoints, 
the switching time required for the magnetization re- 
versal has become of utmost significance today. The 
tremendous increase in storage density and read-write 
speed in magnetic storage media is reaching its limits 
and the need for further speedup of the magnetization 
dynamics is one of the main issues. Actual read-write 
processes operate on the order of nanoseconds and newly 
developed methods such as precessional switching^ have 
demonstrated the possibihty to reach sub-nanosecond 
time scales. Recent pump-probe experiments have also 
shown a fast decay of the magneto-optical signal occur- 
ring on the sub-picosecond time scale using time-resolved 
techniques such as, to cite a few, the magneto-optical 
Kerr effect (MOKE),^ time- resolved second harmonic 
generation (SHG),"^ and pulsed inductive microwave mag- 
netometer (PIMM) 4 

These techniques also make it possible to prepare the 
system in a non-equilibrium state with the magnetization 
vector pointing in a freely chosen arbitrary direction. In 
particular, unlike linear FMR experiments, large-angle 
motion of the magnetization can be studied. This is rele- 
vant to technological applications since magnetic record- 
ing heads are expected to change magnetization direc- 
tions over large angles at ever increasing frequencies. An- 
other advantage of these techniques is that they are sen- 
sitive to spin- wave and damping phenomena. ^iSji 

In this paper we will study the relaxational switch- 
ing of the magnetization initially put in a state far from 
equilibrium with the methods mentioned above. 

In general, relaxation is accompanied by a change of 
the energy of the relaxing system, and the energy re- 
leased or absorbed has to be accommodated by a larger 
system serving as a heat bath. In particular, magne- 
tization reversal of (quasi-single-domain) magnetic par- 
ticles towards equilibrium reduce their Zeeman and/or 
anisotropy energy, and the role of heat bath can be played 
by phonons and/or conduction electrons, etc. Most of 



the existing theories consider magnetic particles as single 
macroscopic magnetic moments whose dynamics is de- 
scribed by the Landau-Lifshitz equation with damping^ 
and Langevin noise field (see, e.g., Refs. S JQJliV 

As was pointed out by Suhli^ and demonstrated by Sa- 
fonov and Bertram,J^i^ internal spin-wave (SW) modes 
in the particle can serve as a heat bath and thus be 
responsible for the particle's relaxation far from equi- 
librium and thereby for magnetization switching. As 
a rule, this internal relaxation process should be much 
faster than those processes that operate via nonmagnetic 
degrees of freedom, with important implications in en- 
gineering of magnetic elements in electronics. Progress 
in computing allowed to simulate magnetic particle's dy- 
namics as that of a system of many interacting spins^'^ 
considered classically. Simulations show that relaxation 
via internal SW indeed occurs. 

From the theoretical viewpoint, SW dynamics in a 
magnetic particle whose global magnetization strongly 
depends on time in the course of reversal is a new and 
challenging subject. Dynamics involving moderate devi- 
ations of the particle's magnetization from equilibrium, 
including the Suhl SW instabilities^^ that may lead to 
a much faster relaxation j^d^d^ can still be described 
with the help of a nonlinear SW theory built around the 
ground state. To the contrast, large deviations require 
redefinition of the spin-wave vacuum by considering SW 
dynamics in the frame related to the instantaneous par- 
ticle's global magnetization. The vacuum corresponding 
to a coUinear state with all particle's spins oriented in 
an arbitrary direction is an example of a false vacuum 
that is unstable and decays toward the true vacuum. A 
spin-wave vacuum of this kind was suggested for the ther- 
modynamics of low-dimensional magnetic systems in zero 
field that have a zero order parameter at T > Its ex- 
tension for nonzero field was done in Ref. [H. Recently, 
the initial stages of SW dynamics including instabilities 
have been considered in the global-magnetization frame 
in Ref. [H. 

The false-vacuum initial condition can be created by 
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biasing the system with a magnetic field that leads to 
the disappearance of the metastable energy minimum in 
which the magnetization is set.^- Other methods are the 
precessional22 and current-induced^^ switching that al- 
low to rotate the magnetization into an arbitrary direc- 
tion during the time that is shorter than the relaxation 
time. 

The aim of this paper, which is the extended version 
of the preceding letter—, is two-fold: 

1. We show by numerical simulations for atomic spins 
on the lattice at initially T — that, for particles 
large enough, excitation of internal spin waves leads 
to a full magnetization switching and relaxation to 
a thermal state with the temperature T defined 
from the energy balance. This solves the puzzle 
of "incomplete relaxation" observed in Ref. [l3| for 
a particle of only 64 effective spins. We show that 
magnetization switching via internal spin waves is 
typically accompanied by a strong reduction of the 
particle's magnetization m that subsequently re- 
covers to a value that is slightly less than the initial 
value m = 1 for coUinear spins. This longitudi- 
nal relaxation is similar to that described by the 
Landau-Lifshitz-Bloch (LLB) equation (see, e.g., 
Refs. HollTll ). although it has a different origin. 

2. In addition to the exponential instability mecha- 
nism that occurs in the case of elliptic magneti- 
zation precessioni^ii^ii^ii^ii^, there occurs a linear 
instability mechanism for the circular magnetiza- 
tion precession. The latter may be driven by a 
random core (volume) anisotropy or by a surface 
anisotropy that causes spin non-coUinearities. We 
analytically derive the rates describing the creation 
of SWs out of the false vacuum, using the spin- 
wave theory in the frame of the particle's global 
magnetization m. Essential dependence of these 
rates on the angle between m and the magnetic 
field H, which shapes the energy landscape, leads 
to a slow non-exponential relaxation of m at large 
times. Furthermore, depending on the direction of 
the applied field, the magnetization switching may 
or may not occur. 

The remainder of the paper is organized as follows. Sec. 
im introduces the classical-spin model of a magnetic parti- 
cle including a bulk anisotropy and a random anisotropy. 
Here the existence of two types of SW instabilities is 
demonstrated. Sec. IIIII contains the formalism of spin- 
wave theory in the frame related to particle's global mag- 
netization. The details of the full analytical solution in 
the case of random anisotropy are given in the Appendix. 
In Sec. llVl the results of numerical atomistic simulations 
are presented that show the two instability scenarios lead- 
ing to the magnetization reversal. Sec. |V] contains our 
discussion. 



II. THE HAMILTONIAN AND SPIN- WAVE 
INSTABILITIES 

A. The Hamiltonian 

We consider the classical Hamiltonian (|si| = 1) on the 
lattice 

T^^^W^j-h-^Sj- JySj-Sj, (1) 

i i ij 

where h — /xqH, /iq is the magnetic moment associated 
with the spin, H is the magnetic field, is the exchange 
interaction, and Ti.Ai is the crystal-field energy at site z, 
a function of satisfying the symmetry of the problem. 
In applications below, we will consider the bulk uniaxial 
anisotropy with easy axis 

nA, = -D{e,-Si)\ D>0 (2) 

and the random anisotropy 

'Hai — - y^^9i,af3SiaSip (3) 
a/3 

with 

gi,al3 = Dr (^UiaUi/3 - ^Saf3^ , Db > (4) 

Ui being a unit vector assuming random directions. — 
One can also add a surface anisotropy and dipolc-dipole 
interaction. Throughout the paper we do not include any 
coupling to the environment. 

The particle's magnetization is defined as 

i 

where Af is the total number of spins. 
Atomic spins obey the Larmor equation 

= [si X fli] , hfli = -dn/dsi. (6) 



B. Spin- wave instabilities 

We now study two models, one with uniaxial 
anisotropy with the same easy axis for all spins and no 
random anisotropy {Db = 0), and the other with ran- 
dom anisotropy and no uniaxial anisotropy (D = 0). We 
will show that in the case of uniaxial anisotropy the SW 
instabilities are exponential while in the case of random 
anisotropy these instabilities are linear. The calculations 
are both analytical and numerical in the latter case and 
only numerical in the former. 
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1. Uniaxial anisotropy and exponential instabilities 

The first model we study is that of uniaxial anisotropy 
with the common easy axis in the z direction, non- 
collinear with the applied field h. In this case, the spin 
precession is non-circular. The same efi^ect may also be 
caused by a biaxial anisotropy. We choose the initial 
state Si = Bx and h =hex- Linearization around this state 
yields the SW spectrum 

Ek = ^J{h + Jo- Jk) {h-2D + Jo- Jk), (7) 

where Jk is the Fourier coefficient of which, for a 
particle with simple cubic (sc) structure, reads 

Jk = 2J^cos(a/cQ,), (8) 

a 

where a is the lattice spacing. In the long-wave-length 
limit Jo — Jk =^ J (a/c)^ . 

Let us now consider the properties of Eq. ((7|) in differ- 
ent regions of the magnetic field h. 

• For h > 2D the state — is the energy min- 
imum and thus Ek is real. This means that there 
is no SW instabilities and the initially excited SW 
will maintain their amplitudes. 

• In the interval < h < 2D the state = e^; is a 
saddle point, so that m can rotate away from this 
state, similarly to the case studied in Ref. Ts'. SW 
modes in the interval < Jo — Jk < 2D — h are 
unstable since becomes imaginaryj^^ This leads 
to the exponential increase of the deviations from 
the initially nearly coUinear state. The highest in- 
stability increment is realized in the middle of the 
k-interval of instability, i.e., Jo — Jk = D—h/2. This 
means that exponentially growing spin waves with 
this nonzero value of k dominate in the instabil- 
ity process. As a result, the magnetization length 
TO = |m| decreases upon rotation out of the saddle 
point, as was observed in early simulations^^ 

• For h < the state = is the energy maximum 
and thus m performs small-amplitude precession 
around e^^ in the absence of SW processes. Indeed, 
according to Eq. ([7]), the mode k = is stable as 
eo is real. On the other hand, k 7^ modes in the 
interval 

-h < Jo - A < 2D - h (9) 

are unstable. In this case the only way of reversal 
is via excitation of internal spin waves with k ^ 
that strongly reduce the magnetization magnitude 
TO. This means that the uniform-precession mode 
decays in favor of non-uniform spin- wave modes. 



S. Random anisotropy and linear instabilities 

The second situation we consider here is with D = 
and non-zero random anisotropy of Eq. ([3]). The latter 
does not break the global particle's magnetic isotropy, 
and its only role is to provide strength for SW conversion 
processes that lead to non-conservation of m-h and thus 
to reversal. We choose h =hez with ft, > 0. Here the 
precession of m is circular and there are no exponential 
SW instabilities. Instead, spin waves may be generated 
out of a false vacuum by linear transformation processes. 
In particular, for m antiparallel to h, spin waves in the 
particle have a negative gap —h. Thus a SW with k 
can be created out of the false vacuum if its energy is 
zero: 

£k = -/i + Jo- Jk = 0. (10) 

In Sec. mil and the Appendix this process will be con- 
sidered in detail for any angle between m and h. The 
amplitudes of unstable SW and the deviation of to from 
saturation increase linearly with time at small times. 

3. Size effects 

In the analysis of both exponential and linear SW in- 
stabilities one has to take into account the fact that be- 
cause of the finite size of magnetic particles their internal 
SW modes are discrete. In particular, for a box-shaped 
particle with free boundary conditions (fbc) there are 
standing spin waves with wave vectors^ 

fca = — — , na = 0,l,...,Na-l, a = x,y,z, (11 
aNa 

where a is the lattice spacing and N^NyN^ = M is the to- 
tal number of spins. For most of other shapes, SW modes 
in magnetic particles have to be found numerically and 
they are labeled by discrete wave numbers rather than 
by the wave vector k. In the sequel, in the analytical cal- 
culations, we will consider only the box-shaped particles 
for simplicity. 

One can see that for the linear instability process, Eq. 
(fTO)) may be satisfied for a particular SW mode. For 
particles small enough, the lowest value of Jq — Jk for 
k 7^ [i.e., for tIq = 1 in Eq. PT|) ] exceeds h and thus 
SWs cannot be created. This yields the absolute stability 
criterion for the particle's linear size L 

L<L* ^ aN* ^ TTflVJ/ft, (12) 

assuming a cubic shape. If there is only one mode that 
exactly or approximately satisfies Eq. (|10[1 and this mode 
does not decay into second-generation SW, there are har- 
monic oscillations between the false-vacuum state and 
the state with the resonant SW mode: 1 — to ~ cos(J7i), 
where depends on the strength of the SW conversion 
processes, i.e., on the random anisotropy or other interac- 
tions creating spin non-coUinearity. These oscillations are 
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similar to the probabihty oscillations between two reso- 
nant states in quantum mechanics. If the unstable spin- 
wave mode can be converted into a second-generation 
spin-wave mode that also has its energy close to zero, the 
relaxation process becomes a two-step process, and the 
temporal behavior of m becomes more complicated. For 
small enough particles, but with L > L* , there are not 
enough resonant SW modes, so that the relaxation pro- 
cess gets stuck in what can be called "spin-wave bottle- 
neck" and the relaxation is incomplete, as was observed 
in ReflH. In large particles, L ^ L*, SW modes become 
quasi-continuous, the SW bottleneck disappears, and a 
cascade of SW processes leads to a nearly full magneti- 
zation reversal, as is demonstrated below. 

Analysis of the exponential instability in magnetic par- 
ticles of finite size goes along similar lines. The particle 
is stable if the smallest value of Jq — Jk with k 7^ (i.e., 
Ua = 1) exceeds the right boundary of the instability in- 
terval, i.e., 2D— /i, see Eq. (jni). This leads to the stabihty 
criterion 



L <L* = aN* 



J 



2D+ \h\ 



for the particle's size L = ma.x{Lx, Ly, L^). Eq. is 
similar to the single-domain criterion since its right-hand 
side is the domain-wall width, if /i = 0. If L > L* and 
there is only one SW mode inside the instability interval 
that does not convert into second-generation spin waves, 
its amplitude and thus the magnetization length m de- 
pends periodically on time. This dependence is not si- 
nusoidal, as that for the linear instability, see comment 
below Eq. mi). Instead, the SW amplitude initially 
exponentially increases but then this evolution becomes 
inverted due to nonlinear effects at large amplitudes, and 
the SW amplitude reversibly returns to the starting small 
value. After that the process repeats periodically. If 
second-generation unstable spin waves are created, the 
time evolution of m becomes more complicated and m 
does not return to saturation. Still, for L > L* the mag- 
netization reversal is incomplete, as observed in Ref. Il3l . 
Only for large enough particles, L ^ L* , does relaxation 
via internal spin waves lead to a complete magnetization 
reversal. 

Finally, we address the question as to how large the 
particle must be so that its spectrum of spin waves be- 
comes effectively continuous. In the case of the exponen- 
tial instability the criterion is that there must be several 
modes within the instability interval whose length (in en- 
ergy) is proportional to D, as we have seen above. Thus 
the criterion of the quasi-continuous spectrum becomes 



Ae < D, 



(14) 



where Ae is the average distance between the SW modes. 
In the case of random anisotropy, the spin-wave devia- 
tions, associated with a particular mode, satisfy linear 
differential equations with a source that is proportional 
to the random anisotropy constant Dn, see e.g., Eq. 



(IA23p . The response to the off-resonance source is of 
order Dji/Ae that becomes large for Ae < Dfi, formally 
coinciding with Eq. p^ . If this condition is fulfilled, 
spin waves are effectively generated for all values of the 
field h without the necessity to satisfy Eq. pO)) . Next, 
Ae can be expressed via the density of SW states, defined 
by Eq. ((J25|) . as 



Ae = 



(15) 



in the vicinity of the energy e = huj. Thus Eq. (|14p can 
be rewritten in terms of the particle's size as 



(16) 



Adopting the rough estimation h/ p{lo) ^ J over the 
whole range of SW frequencies, one obtains the criterion 
of the quasi-continuous spectrum in the form 



(17) 



(13) and the same with D Dn for random anisotropy. 



III. ANALYTICAL THEORY OF SPIN- WAVE 
INSTABILITIES IN THE PARTICLE FRAME 

In this section, we present our general formalism of 
SW theory in the frame related with the particle's global 
magnetization for an arbitrary direction of the applied 
field. Then, we study the two anisotropy models of sec- 
tions [TIBT] and [HB^ and the ensuing SW instabilities. 

The microscopic effective field ftri; in Eq. ([6]) can be 
written in the form 



an. 



dsi 



(18) 



where Ti.Ai contains only non-random anisotropy, 
whereas the random anisotropy is singled out. In particu- 
lar, for Hai given by Eq. ^ one has —dHAi / dsi— D Sj = 
< — > < — > 
D • Si and the components of the tensor D read 

(D)q^ = DSaz^pz- Similarly, the components of the 
random-anisotropy tensor are given by (gT)Q/3 = gi,ai3 
and Eq. One can represent in the form 



m 



(19) 



where m is the average spin defined by Eq. (O and 
■01 contains the Fourier components with k ^ and de- 
scribes spin waves in the particle. Whereas in the stan- 
dard SW theory m is a constant corresponding to the 
ground-state orientation, here it is treated as a time- 
dependent variable. Since the atomic spins are sub- 
ject to the chiral constraint = 1, one can usei^ 
m = n-y/1 — ifJj with n. ■ -^Ji = 0, where n is a unit vec- 
tor. Although this reduces to two the number of the -i/'i 
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components to deal with, the formahsm becomes much 
more cumbersome, the final results, however, are not af- 
fected. Thus we decided not to use the chiral constraint 
explicitly in our presentation. Of course, properly writ- 
ten equations must satisfy this constraint that can be 
used to check them. 

The equation of motion for m following from Eq. ([5]) 
has the form 



hm = [m X hcff] -I- R. 



Here 



hcff = h - 



^nA^ 



dsi 



(20) 



(21) 



is the effective field acting on the particle as a whole and 
does not contain the random anisotropy and exchange 
coupling. The term R couples the dynamics of m to 
that of spin waves described by ipi and it is responsible 
for the relaxation of m. Calculation yields 



R 



-y 



{[mx2'g'iVi] 



[t/.. X ( r'V'. + 2't^ (m + Vj)] } • (22) 



Here [mx2^im] = as an average of the random 
anisotropy, ^j^- [-0^ x V'j] = by symmetry, while 
some other terms vanish by virtue of Eq. (|19p. The 
tensor F is given by 



a/3 



(23) 



For Uai given by Eq. ^ one has [¥] = 2(H)a/3 = 

V / a/3 

In turn, the equation of motion for can be obtained 

as 



ii^ = Si-m = [si X Jl.,] -m. 
Working out the various terms yields 



hij^i = [mx2'g'im] + Af ^ + A 



,(2) 



(24) 



(25) 



where A^^-* contains terms linear in i/jj and not contain- 
ing *g*j, while contains the terms of order t/i^ and 

■0g. In the sequel we will keep only A.^^-* that is respon- 
sible for the generation of spin waves out of a false vac- 
uum, whereas AJ responsible for nonlinear spin-wave 
processes will be dropped. One has 

A^ = X (heff + Jom)] 



mx F 1/"^ + ^ JijTpj 



(26) 



The first term in Eq. ([25)1 that causes non-coUinearity 
of spins, is responsible for the linear SW instability. The 
same effect is produced by surface anisotropy and dipole- 
dipole interaction. This source term in the linear equa- 
tion for "i/jj is due to the lack of translational invariance 
of the random anisotropy. On the contrary, the terms 
due to the bulk anisotropy TiAi and magnetic field h are 
invariant by translation and thus enter the equation of 
motion for the particle's global magnetization m ((20|) 
rather than the equation for -0^. 

In Eq. ([20| it is convenient to project the relaxation 
term R onto m and the perpendicular directions. For 
this we introduce orthogonal unit vectors 



ei 



m 

m 

[n X hcff] 



[n X hcff] 



62 = [n X ei] = 



n X hcff] I /loffVl - x'^ 
[n X [n X heff]] 



fn X h, 



off 



where 



n • h. 



■cff 



^cff 



mhcs 



R then has the form 



R =i?lin-|-i?iei -t-i?2e2. 



(27) 



(28) 



(29) 



where describes the longitudinal relaxation of the par- 
ticle's magnetization (change in the magnetization mag- 
nitude) while 



R2 = Ra 



(30) 



describes the transverse relaxation (rotational relaxation 
of the magnetization vector), since 62 • n = 0. On the 
contrary, Ri does not describe any relaxation; it merely 
describes a small modification of the particle's precession 
due to excitation of spin waves, an effect that will be 
neglected here. Explicitly Eq. ([20]) can now be rewritten 
as 

m [m X [m X hpffll , , 

;im= mx heff ^ r-^, 31 

" m m m X heft 



where 



Rh — n • R, 



(32) 



Eq. ((3T|) resembles the Landau-Lifshitz-Bloch (LLB) 
equation since it includes both the longitudinal and 

f2l 

transverse relaxation terms. Ignoring A^- in Eq. (j25p . 
one can solve the resulting linear equation for tpi, insert 
the solution into R, and obtain the relaxation terms in 
Eq. ((5T|) from Eq. ([5^ . Integrating out V'i can be done 
in the particle's frame defined by Eq. (|?7|) . It is under- 
stood that Eq. ((3T|) is only valid during the initial stage 
of the evolution out of the completely or nearly coUinear 
state, |m| = 1, when the deviations 4'i due to spin waves 
are small and can be considered perturbatively. 
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Linear instabilities 



Eq. ((25)) without the contribution can be solved 
analytically in the case of a pure random anisotropy, 
'^Ai = 0, since then the time dependence of the frame 
vectors is a simple precession around the external field 
h. The details can be found in the Appendix. For I? = 
the dependence in{t) [neglecting R in Eq. ([20]) ] is a cir- 
cular precession, and thus the calculations simplify and 
lead to the equation of motion 

1 r ,1 / N m ■^/9 , , fmx [m X hll 

m = - [m X h]-m^^^Tn{x) m^'^T^{x)- ^— -. 

n m m'^h 

(33) 

Here the transverse and longitudinal relaxation rates de- 
pend on the orientation of the particle's magnetization 
vector X = (m • h) / (mh): 



Til (a;) 



Vf'^„(^)-r„o$„(.) (34) 



TA^) = ^^^<^>±{x)^r^o^^{x) (35) 



with 



$11 (a;) = 



{l~x)' 



4 

1 - x 
6 



{1 + 2xf + V2 {1 - x'^)] (36) 



(1 + 2xf + %/2(l - x)(2 + x)] (37) 



The applicability of our method requires T±, Fn <^ ujh = 
h/h, i.e., the relaxation of the magnetization (as well 
as the rate of SW production) is much slower than the 
magnetization precession considered as unperturbed in 
the first approximation above. 

We note that Eq. ([55]) is similar to the LLB 
equationii^iii However, it is valid, in general, only for 
short times, when the number of excited spin waves is 

(2) 

still small and m = 1, so that neglecting A^ ' in Eq. 
(1^ is justified. Correspondingly, m in Eq. may be 
replaced by 1. Of more consequence, however, is to intro- 
duce the magnetization direction n and write m — nm. 
Then Eq. (|33p can be split into the two following equa- 
tions 



1 r , -, / ^ fnx [n X hll 
-[nxh]-F^(a;)l L_ 



where x = n ■ h/h, and 



The small-t behavior 



-F||(x). 



m,{t) = l-r||(-l)i 



(38) 



(39) 



(40) 



following from Eq. ((39l) and shown in Fig. |4^. It agrees 
well with the numerical result with a small discrepancy 
stemming from h/J not being small enough to use the 
analytical expression for the density of states given by 



Eq. (jA27p . Equation. ([39)1 does not describe the increase 
of m after switching and its recovery to to = 1 that is 
seen in Fig. U^. 

A striking feature of our result is that both F^(a:;) and 
r|[(j:) vanish for m||h (i.e., for x — 1), while they reach 
their maxima at x — —1. Thus, initially fast relaxation 
slows down when the particle approaches equilibrium, see 
Fig. m Indeed, Eq. ([55)) can be rewritten as 



(41) 



x = Tjl{x){1-x^). 



In terms of the angular deviation from equilibrium we 
have 



y ■ 



1 



■ h./h = 1 - x < 1 



(42) 



upon which Eq. (|4ip simplifies into the equation ij = 
— 3r^o2/^ and whose solution reads 



y(o) 



l + 3y(0)rioi' 



(43) 



It is seen that the long-time asymptote of this solu- 
tion does not depend on the initial condition j/(0) and is 
given by y{t) — 1/ (SFj^oi) . The full magnetization vec- 
tor m ~ nm includes m that follows from Eq. (|55)) . The 
latter becomes to = -(9/4)F||o?;^ = -(3/2)F^oy^ = 2;/2. 
This yields 



1 - m(t) = - [y(0) - y{t)] 



(44) 



Thus, we see that the deviation 1— m{t) remains fi- 
nite and small for t ^ oo since y{t) — > 1/ (3r^o^), that 
is the consequence of the relaxation slowing down as m 
approaches equilibrium. The change of m due to its rota- 
tion and due to the change in its magnitude are compa- 
rable with each other. In addition, one has to remem- 
ber that the equation of motion for m was obtained 
for the initial stage of relaxation only and, in partic- 
ular, Eq. (|55)) does not describe thermalization. Still 
the evidence provided by these analytical calculations of 
the non-exponential relaxation via internal spin waves in 
magnetic particles is quite convincing. 



B. Exponential instabilities 

In the case of the purely bulk anisotropy (no random 
anisotropy and other interactions causing noncoUinearity 
of spins), Eq. (|^5)) with A^^^* is a system of uniform 
linear equations. In general, m has a nontrivial quasi- 
periodic time dependence defined by the bulk anisotropy, 
magnetic field, and the energy associated with m that is 
conserved in the absence of the internal spin waves. Ex- 
citation of the latter reduces the energy associated with 
m, so that the total energy is conserved. Mathemati- 
cally spin-wave instabilities in this case are exponential 
divergences of the solution of a system of linear differ- 
ential equations with periodic coefficients. A well-known 
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example of such instabilities is parametric resonance. In 
all cases instability requires elliptic or more complicated 
precession of m. Simple circular precession around the 
field h does not lead to exponential instabilities. 

From the physical point of view, exponential SW in- 
stabilities are similar to the Suhl exponential instabili- 
ties. The difference resides mainly in the mathematical 
description. The Suhl formalism considers spin waves 
above the true ground state, so that it cannot be applied 
in situations of large deviations of m from the energy 
minimum. The advantage of this method is that unstable 
spin waves are described by a linear differential equation 
with constant coefficients that can easily be solved. 

In contrast, the present method using the frame related 
to the particle's global magnetization can be used for any 
deviations of m from the ground state. The price to pay 
is to deal with differential equations with time-dependent 
coefficients that are difficult to solve analytically (see also 
discussion in Ref. .19,) . However, if m is initially oriented 
towards an energy maximum or a saddle point so that it 
does not evolve in time in the absence of spin waves, one 
obtains differential equations for the deviations ipi with 
constant coefficients that can easily be solved showing 
exponential instabilities. In fact, in this case one can just 
linearize Eq. ^ near a given direction of m, instead of 
going through the formalism of this section. Examples 
are considered below Eq. ([7|). 



IV. NUMERICAL RESULTS 

Efficient magnetization reversal via the SW instabil- 
ity requires a quasi-continuous SW spectrum, i.e., the 
existence of many SW modes in the instability interval. 
For this, a strong inequality opposite to Eq. (jl3[) should 
be fuffiUed, i.e., the particle's size and/or the anisotropy 
D should be large enough. To illustrate the process, we 
have solved Eq. ^ for magnetization switching using the 
C-I--I— based package MagneticP article (© H. Kachkachi 
& L. Reynaud). 

For the model with uniaxial anisotropy with easy axis 
along z, the calculations have been done for two differ- 
ent particle's shapes. One is a sphere of radius i? = 35 
lattice spacings, cut inside a 70 x 70 x 70 cube resulting 
in 171712 atomic spins. The other particle is a paral- 
lelepiped consisting of 28 x 31 x 34 = 29512 atomic spins. 
The idea behind choosing all different sizes is to avoid de- 
generacy in the energies of SW modes [see Eq. (ITT]) ] and 
ensure a smoother density of states and thus facilitate 
SW conversion processes. In both cases the initial state 
is coUinear and a very small surface anisotropy was added 
to create initial very small deviations from the coUinear- 
ity that then exponentially grow. The lattice structure 
is simple cubic and boundary conditions are free (fbc). 
The magnetic field is h = he^, with h > 2D that creates 
the energy maximum in the direction opposite to h. The 
latter is the initial direction of the magnetization in part 
of the simulations. As stressed above, a particle whose 




Figure 1: Magnetization switching out of the maximal-energy 
state via exponential spin-wave instability in a spherical par- 
ticle with uniaxial anisotropy and transverse field. Switching 
occurs via longitudinal relaxation without rotation, (b) shows 
the short-time magnification. The "Plateau" at t < 2000 de- 
scribes the initial exponential increase of 1 — m. 

magnetization is oriented toward the energy maximum 
cannot just rotate out of this orientation (in the absence 
of a coupling to the environment) because of the energy 
conservation. 

The results shown in Fig. [1^ for the spherical parti- 
cle and in Fig. for the box-shape particle are similar. 
One can see a nearly full reversal via the longitudinal re- 
laxation, since the transverse magnetization component 
mj_ remains small at all times. In the middle of the 
switching process, excited spin waves almost completely 
destroy the magnetization m. Then the first-generation 
long- wavelength SWs of high amplitude convert via non- 
linear processes into all possible spin wave modes and 
the system thermalizes. Since the thermodynamics of 
classical spins is mainly determined by short-wavelength 
modes of high energy, the energy conservation requires 
that the amplitudes of these SWs be small. This ex- 
plains the almost full recovery of m after reversal. The 
asymptotic disordering 1 — rn > corresponds to the final 
temperature T following from the energy balance. 

The short-time magnification in Fig. shows an ex- 
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Figure 2: Magnetization switching via exponential spin-wave 
instability in a box-shaped particle with uniaxial anisotropy 
and transverse field, a) The particle is prepared with all spins 
opposite to the magnetic field (the maximal-energy state). 
Again the longitudinal relaxation is the main mechanism of 
switching, b) The particle is prepared with all spins perpen- 
dicular to the magnetic field. 

ponential increase of 1 — m due to the SW instability until 
the first dip, followed by an incomplete recovery of m and 
its further decrease. The same behavior was observed for 
other particle's shapes. 

If the initial particle's magnetization makes an angle 
with the highest-energy direction, the particle's dynamics 
is a combination of a non-circular precession of m and re- 
laxation via SW processes, so that our numerical results 
show a more complicated behavior. The closer is the ini- 
tial orientation of m to the energy minimum, the slower is 
the relaxation towards it and the smaller is the decrease 
of m in the course of relaxation. The results for the 
box-shaped particle with its initial magnetization vector 
perpendicular to the highest-energy direction are shown 
in Fig. ^Bp. One can see that the transverse relaxation is 
extremely slow in this case, so that a prohibitively long 
computer time is needed to follow the decay |m^| —t 0. 

Fig. [3] shows the magnetization reversal out of a dis- 
appearing metastable state for the model with uniaxial 
z anisotropy and oblique magnetic field applied at the 
angle -0 = ''^/^ to the z axis. For h = D there is a 




Figure 3: Magnetization switching out of the disappearing 
metastable state via exponential spin-wave instability in a 
particle of a box shape with uniaxial anisotropy and oblique 
magnetic field. 



disappearing metastable minimum for the orientation of 
the spins at 6* = 37r/4 to the z axis. Indeed, the energy 
E{0) = -Dcos^e- {h/V2) cos 9- {h/V2) sin 61, follow- 
ing from Eqs. m and 1^, satisfies E'{0) = E"{e) = for 
h = D and 9 = 37r/4. Also one can check that for h = D 
one has E'{9) = at 61 = 7r/12 that is the energy mini- 
mum. Notice that for simulations, a slightly higher value 
of h has been chosen to let m precess away from the ini- 
tial state with 9 = 37r/4. One can see that this precession 
is strongly non-circular and the particle's global magneti- 
zation m tends to return into the initial state that would 
comply with the energy conservation law for a single spin. 
However, due to the excitation of internal spin waves in 
the particle the energy associated with m decreases, so 
that the returning to the initial orientation is incomplete. 
After a couple of cycles a lot of exponentially unstable 
spin waves get excited and a faster relaxation of m be- 
gins. Asymptotically m approaches the orientation close 
to the energy minimum at 9 = n/12, up to thermal dis- 
ordering. The reduction of m increases exponentially at 
small times but does not exceed 40% in the middle of 
relaxation in this case. 

Fig. |4] shows the results of our numerical simula- 
tions for the magnetization switching in the random- 
anisotropy model described by Eq. ([3]). The particle is 
again of a box shape with fbc, consisting of 28 x 31 x 34 = 
29512 atomic spins. One can see again a nearly full rever- 
sal which is, in contrast to Fig. [TJ), linear at short times. 
The simulation results at short times are in a reason- 
able accord with the analytical result in Eq. that 
is shown by the dashed line. Also shown in Fig. [4]i is 
the asymptotic approach of the thermal state with tem- 
perature T corresponding to the released Zeeman energy. 
Note that the reduction of m in the middle of reversal is 
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Figure 4: Magnetization switching via linear spin-wave insta- 
bility in a particle with random anisotropy. a) Particle pre- 
pared with all spins antiparallel to the magnetic field. The 
small-t asymptote of Eq. (|40|l is shown by the dashed line. 
For this initial condition, switching occurs predominantly via 
changing the magnetization length m. b) Particle prepared 
with all spins perpendicular to the magnetic field. 

much stronger for the anticoUinear initial state (Fig. 0^) 
than for the perpendicular initial state (Fig. \^), simi- 
larly to the case of the exponential instability. 

As discussed in Sec. [TTl in smaller particles SW modes 
are essentially discrete and SW conversion processes are 
blocked by the impossibility to satisfy the energy con- 
servation law. Even if one SW mode is generated as a 
result of exponential or linear instability, it usually can- 
not be converted into other modes. This leads to the 
spin-wave bottleneck in the relaxation process via inter- 
nal spin waves. As an example we consider the sc-lattice 
cubic particle with = Ny = = h [N = 125) the 
three degenerate modes (1, 0, 0), (0, 1, 0) and (0, 0, 1) are 
excited via the linear instability process in the model 
with random anisotropy for the magnetic fields in the 
vicinity of h/J = 0.38197, as follows from Eqs. dH]), 
and pT|) . In Fig. [5] obtained with Dn/ J = 0.01 one can 
see that at h/J — 0.38197 there is an essential reduction 
of the initially saturated magnetization, although the re- 
laxation is bottlenecked because second-generation spin 




Figure 5: Time dependence of the magnetization magnitude 
m for the particle of a cubic shape consisting of 5 x 5 x 5 = 125 
spins with random anisotropy. For the resonant value of the 
field h spin waves are generated out of the false vacuum (mag- 
netization opposite to the field) but they cannot be converted 
into other modes so that the process is bottlenecked. For the 
non-resonant value of h spin waves are practically not gener- 
ated. 

waves cannot be created. The evolution of m is not si- 
nusoidal because three degenerate SW modes are excited 
with different efhciencies, depending on the realization 
of the random anisotropy, so that each mode has its own 
conversion frequency fi, defined below Eq. p^ . To the 
contrast, for the off-resonance value hj J — 0.37 there 
is no SW mode at resonance with the false vacuum and 
deviations 1 — to remain very small. 

For larger values of the random anisotropy the par- 
ticle's SW spectrum becomes effectively continuous, in 
accordance with the arguments at the end of Sec. |TT1 In 
particular, for D^/ J = 0.1 there is no principal difference 
between the results for h/J = 0.38197 and h/J — 0.37 
and the magnetization relaxes towards the ground state. 
However, this relaxation is much slower than in really 
large particles considered above and there is a significant 
quasi-random dependence of the relaxation rate on model 
parameters. 



V. DISCUSSION 

We have shown that internal spin-wave processes in 
large enough magnetic particles can lead to complete 
magnetization reversal and thermalization. The role of 
thermal bath in these processes is played by the magnetic 
particle itself. The energy release in the course of the re- 
laxation of the particle's global magnetization m towards 
the ground state is absorbed by internal spin waves, so 
that the total energy is conserved, as long as there is no 
coupling to the environment. 
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The two main scenarios of magnetization reversal are 
through exponential and linear SW instabilities. While 
the former is the Suhl instability, described here within a 
more general formalism that allows for large magnetiza- 
tion motions, the linear instability is difficult to pinpoint 
in the existing literature. Anyway, theoretical descrip- 
tion of the linear instability requires redefinition of the 
spin-wave vacuum by using the frame related to the par- 
ticle's global magnetization m that in general depends 
on time. 

In both cases, the relaxation is fast at the beginning 
but then slows down. These results are relevant to the 
study of the dynamics of magneto-electronic elements 
and, more generally, in the physics of unstable macro- 
scopic states. In the model with random anisotropy 
that can be solved analytically relaxation at asymptot- 
ically large times is power law rather than exponential. 
The reason is that the relaxation rate becomes small as 
the direction of the particle's magnetization approaches 
the ground-state direction. For the random-anisotropy 
model both longitudinal and transverse relaxation rates 
are maximal for m pointing toward the energy maximum. 
It is clear that the majority of the magnetic nanoparticles 
studied nowadays are shown to exhibit effective uniaxial 
anisotropy. However, as their size is further reduced, 
the on-site crystalline anisotropy may adopt a more dis- 
ordered distribution upon which the model of random 
anisotropy may become of some relevance. 

Estimations of relaxation rates via internal spin waves 
in metallic Co made in Ref. [H yield T -10^ - 10*^ a~'^. 
Comparison with the spin-phonon and other external 
rates, obtained microscopically, is difficult since the lat- 
ter are still unreliable and should be modified by collec- 
tive processes such as phonon bottlenec k^^'^^ and phonon 
superradiancei^ In most cases internal processes should 
dominate. One should not forget, however, that the lat- 
ter die out for nanoparticles that cannot accommodate 
spin waves. The phenomenological Landau-Lifshitz re- 
laxation rate T^l = ah/h with the damping constant a 
ranging between 10"^ and 10^'^ is typically much larger 
than microscopic rates that do not take into account col- 
lective processes. One has to be careful, however, since 
the a that is extracted from experiments may contain a 
contribution from the SW processes. 

Numerical simulations in this work have been done for 
a classical-spin model without coupling to the environ- 
ment. The dynamical equations are Larmor equations 
for each atomic spin precessing in the effective field cre- 
ated by other spins plus the magnetic and anisotropy 
field. Including the effect of the environment by a dou- 
ble vector product relaxation term introduced by Landau 
and Lifshitz^ is straightforward. Also one can include the 
finite-temperature effect by the Langevin fields that will 
slow down the program, however. It should be stressed 
that such a method of including the environment is not 
reliable because it misses collective effects in the spin- 
lattice relaxation. 
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Appendix A: RANDOM ANISOTROPY: 
DYNAMICS IN THE PRECESSING FRAME 



Consider a model with Ti-Ai — and small random 
anisotropy *g*i that is the only source for the particle's 
relaxation. In the zeroth approximation m is simply pre- 
cessing around the magnetic field, 



as well as 



hm = [m X h] , 



(Al) 



and it is convenient to consider spin waves in the frame 
related to m and defined by Eq. ((27l) . Eqs. ((27|) can be 
solved for h resulting in 



(A2) 



h =hxn—h\J 1 — x'^e2 



Time derivatives of the basis vectors defined by Eq. (127] 
are given by (using (jAip ) 

hh — [n X h] 

[hh X h] [[n X h] X h] —nh + hx 



hei = 



he-2 = 



hh (n • h) [n X h] : 

hy/l — x' 



(A3) 



and they can be projected on n, ei, and 62 : 



hh — h\/ 1 — x^Gi 

hei — (n • hei) n + (e2 • hei) 62 = —h\/ 1 — x'^n — hxe2 
he2 — hxei. (A4) 

Now the explicit time dependence have to be obtained. 
To this end, it is convenient to choose the z axis along h, 
i.e., h = hez, then from Eq. (jAl|) follows 



n{t) = xe^ — — x^ cos{uJht)ex + \/\ — x^ sm{ujht)ey 

ei{t) = cos{uJht)ey + siii{LOht)ex 

e2{t) = — Vl — x'^Bz—x cos{uJht)ex + xsiii{uJht)ey. (A5) 
where 



ujh = h/ h. 



(A6) 



One can project tpi onto the time dependent n, ei and 
62 as 



ipi=ipMn + ipiiei + Tpi2e2, 



(A7) 



and thus 



hipi = hipion+hipiiei + hipi2e2 

+ Tptohh+il^iihei + ilJi2he2. (A8) 

(2) 



Inserting this form into Eq. ((25)l with TlAi — and A- 
dropped and projecting onto the three basis vectors one 
obtains the equation 



hipio — ipiih^/ 1 — x^ = n • A, 



(1) 



htpii + ^pioh\/l - a;2 + hx^pi2 
= ei • [mx2*g^im] + ei • a|^^ 

together with 

fiipt2 - hxipii = 62 • [mx2^im] + 62 • Af ^ 
In the rhs of Eq. (jA9p one has 

n • A^^ = n ■ [ipi X h] = • [h x n] . 
Then with the help of Eq. ((27l) one obtains 
V'io = 0. 



(AlO) 
(All) 
(A12) 
(A13) 



Expressions in the rhs of other equations can be processed 
as follows 

ei • [mx^im] 
= m^ei • [nx {ei (ei^in) +62 (e2*g'in)}] 
= m^ei • {e2 (ei^jn) - ei (e2'g'.in)} 
= -'m^{e2^^n) (A14) 

and, similarly, 

62 • [mx^jm] = w? (ei^g^^n) . (A15) 

From Eq. ((26l) with T-Lm — one obtains 



Af ^ = X (h + Jom)] + mx ^ J^Vj 

i 

= [(■0ion + V'jiei +^^^262) X (h + Jom)] 

+ mx ^ (V'jon + Vjiei + -0^262) 
3 

= ipio [n X h] + (Jei x h] + Jqto Jei x n]) 
+-012 (le2 X h] + Jqui [62 X n]) 

+m J,j (tpji [n X ei] + tpj2 [n x 62]) (A16) 
j 

With the help of Eqs. ^7}i and ((J2|l this becomes 

Af' = iPmhVl - x^ei 

+"011 (^—hxe2~h\/ 1 — x^n — Jome2^ 
+tpi2 {hxei + Jomei) 

+ra^ Jij (-0^162 - ^pJ2el) (A17) 
i 

whose components of this vector are 

A\1^ = hxtpi2 + (Jo% - Jy) 0j2 



(A9) 



A, 



(1) 

i2 



+ h^/l - x^ipio 
'hxtpa - (Jo(5.y - Jy) V'ji. (A18) 
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Now Eqs. (jAlOp and (|A11[) after simplification become 

hipii = (Jo(5y - Jij)'4'j2 - 2m^ (e2^jn) 

j 

h'ii)i2 = ~m^{JoSij - Jij) Tpji + 2m^ (ei*g*in) . 



This is the system of equations describing the generation 
of internal spin waves in a magnetic particle. 
It is convenient to introduce the variables 



for a particle of a box shape with fbc for huj <^ Jq one 
has 



Pex(w) 



a 



n (WcX,k - 



k dfc(5(ecx,k - e) 



2n f V 



(A27) 



V'»,± = ±iV'i2, e± = ei±ie2 (A19) Eq. (H23| has the solution 



so that 



1 i 
ei = -(e_+e+), 63 = - (e_ - e+) 

1 i 

and 

eiV'ii + e2V'j2 = i?£ (e_V-'i,+) 

ei'(/'i2 - e2V'ji = /r7i(e_'0i,+) 

ei'(/'j2 + e2'(/'ji = (e+V-'i,+ ) • (A20) 

Then one obtains the equation 

h4'i,+ = -zm^ (Jo<5.y - Jij) + 2im^ (e+^in) . 

(A21) 

what about ipi^-l 

For a box-shaped particle one can rewrite this equation 
in terms of the discrete Fourier components 

i k 

etc., as 

fi-0k,+ = -i£cx,k^Ak,+ + 2zm^ (e+'g^kn) , (A23) 
where 

£ox,k = m ( Jo - Jk) . (A24) 
The density of pure-exchange spin-wave states is given 

by 



Pox(t^) = -J^^5 (^cx,k 



UJ] 



(A25) 



where fiucx.k = £cx,k and M — N^NyNz- It satisfies 
/ dujp{uj) = 1. In the continuous approximation for small 
wave vectors for the sc lattice one has 



eox,k = mJ {ak) 



(A26) 



^k,+ W = ^k,+ (0)e 



+ 2m^- I dt' 



xe 



(e+(t')g^kn(t')) 



(A28) 



where ?iWex,k = £ox,k- The first term of the solution takes 
into account spin waves already available in the particle, 
such as thermal spin waves. The second term describes 
spin waves generated by the particle's precession via the 
random anisotropy. Note that, according to Eq. (jA13[) . 
longitudinal fluctuations do not have any dynamics 
in the linear approximation. 

The next step is to substitute the solution for Vk,-i-(0 
into the relaxation term R in the equation for m, Eq. 
(f20| . Keeping terms of order V'g in Eq. ([22]) with TiAi — 
one obtains 

R = 2m-^ { [nx'g',V^»] + [^H X "g^^n] } . (A29) 



For i?|| of Eq. (I32|) . using transformations in Eq. (jAl4p 
one obtains 

i 

■i 

= ~2to-^ ^ -0r [n X (ei (ei'g^in) + 63 (es'g^in))] 

■i 

= -2?7i-^^V»- {e2 (ei'g'^n) - ei (e2'g'jn)} 

i 

= 2m^Y{{e2^irL)'ipii- [ei^ir\)il>a] ■ (A30) 



Finally, with the help of Eq. ((X20l) 



R\\ = —2m Im 



(A31) 
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Further one obtains 

i?2 = e2-'R.=2m^^e2- {[nx'g'iilji] + x'^iu]} 

i 

= 2m-^ ^{{ei'^t^pi) - A ■ [e2 x ^,n\ } 

i 
i 

- Jpi ■ (ei {n^.n) - n (ei'g'.n)) } 
= 2m-^ ^ {(e2'g*iei) -^ii + (e2*g*ie2) V'j2 

- Vzi(n'^^n)}, (A32) 

where we have dropped all the terms with "i/jio since they 
are frozen in and disappear after the averaging over the 
random anisotropy. This can be rewritten as 



i?2 = 2m Re 



r? { (ei'g'^e^) - (n^^.n) } ^,,+ 



(A33) 

In the Fourier representation R\\ and i?2 = -R± are 
given by 



-Rii 



-2m /m 



i?i = 2mi?e 



/2 X! 



1, •/ 



X {(ei'g'^ke-) - (n'g'^kn)} Vk,+] (A34) 
Note the relations 

"g^k 'g'-k, ^k,± = V'-k,T- (A35) 
Using Eq. (|M8)) yields 

X (e_(t)^_knW) (e+(t') "gknlt'))] •(A36) 
Random anisotropy has the form of Eq. @ , so that 

3k,Q/3 = -Di?, ^ e'-^'^'UiaUip. (A37) 

Using the formula for the ensemble averaging on a lattice 
site i 
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— {SafjSjS + Saj6ps + l5a5<5/37) + gSapS^S, (A38) 



whereas the correlator of two g on different lattice sites 
is zero, one obtains the important formula 

-J^((a^kb) (c^-kd)) 



Af 
15 



{(a • b) (c • d) + (a • c) (b • d) + (a • d) (b • c)} 



-^(a-b) (cd) 



(A39) 



In fact, the terms with (a • b) (c • d) vanish in the expres- 
sions below. With the help of Eq. (jA39p one obtains 



n2 
15 



l.yl o.,.(t-0 



x{(e_(t).e+(t'))(nW-n(0) 
+ (n(t).e+(0)(e_(t)-n(0)}]. 



(A40) 



After computer algebra using Eqs. (jA5p and (|A19p one 
obtains 

(e_(t).e+(t'))(n(t)-n(t')) 
+ (n(t).e+(0)(e-(t)-n(i')) 
= (1 - x^) 

+ i(l-a;)'(l + 2x)'e""" 

+ ^{l-xf{l-x^)e^'M'-'') 

+ 1 (1 + xf (1 - x^) e-2»'^4*-*'), (A41) 



where x is given by Eq. (^5]) and w/j by Eq. (jA6|) . Since 
in Eq. (|A40|) Wox,k > 0, here one should keep only the 

resonant terms with e*"''(*"*') and e^''^'^(*"*') that sat- 
isfy the energy conservation and cause transitions. One 
obtains 



15ft AA^ 



X <{ i (1 - x)^ (1 + 2.x)^ 7r(5 (tJox,k - t^/i) 



+ i (1 - X)' (1 - X^) ^5 (L^e 



-— — ^(l-a;)' (l + 2x)VpexK) 
15n L 



+ (l - x^) 7rpox(2wh)] . 
With the help of Eq. (|A27|) one obtains 



Rn 



157r J V J 
where $11(2:) is given by Eq. ([M)) 



(A42) 



(A43) 
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For i?2 from Eqs. (|A34|) and (|A28p one obtains 



This results in 



X {(ei(i)"i'_ke-W) - (n(i)^_knW)} 
X (e+(t')¥kn(t'))] 



^"1 



15 



1 X - i 



x{(ei(0-e+(t'))(e-(t)-n(O) 
+ (ei(i)-n(0)(e-(t)-e+(0) 
- 2(nW.e+(O)(n(0-n(O)}] 



(A44) 



Computer algebra yields 



{(ei(t).e+(0)(e-(t)-n(i')) 



+ (ei(t).n(t'))(e-W-e+(t')) 
- 2(n(<).e+(0)(n(t)-n(t'))}] 

-3x3 - 1(1 - a;)(l + 2a;)2e*""(*-*') 
+ i(l + a;)(l-2a;)2e-*'^'-(*-*') 
-i(l-x)2(2 + a;)e2*'^''(*-*') 
+ i(l + 2;)2(2-a;)e-2»-"'(*-*'). 



277i3 £)2 

i?± = T^Vl ~ - x)tt {(1 + 2a;)VcxK) 



15/1 

X (l-a;)(2 + a:)pcx(2t^?0}- 



(A46) 



Finally, using Eq. (jA27P again, this expression can be 
brought into the form 



Ra 



v"^Dl [h 



57r J 



^yfyT^$±(x), (A47) 



where is given by Eq. (|37|1. Finally, with 



T,{x) = -\r,{.1 r,(x) = -i-^ (A48) 



(A45) one obtains Eq. ([33 
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